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Abstract 

We make a QCD light-cone sum rule (LCSR) assessment of B( g \ semi-leptonic decays to a light 
scalar meson, B^ — > SWi,Sll (I = e,u,r). Chiral current correlators are used and calculations 
are performed at leading order in a s . Having little knowledge of ingredients of the scalar mesons, 
we confine ourself to the two quark picture for them and work with the two possible scenarios. 
The resulting sum rules for the form factors receive no contributions from the twist-3 distribution 
amplitudes (DA's), in comparison with the calculation of the conventional LCSR approach where 
the twist-3 parts play usually an important role. We specify the range of the squared momentum 
transfer q 2 , in which the operator product expansion (OPE) for the correlators remains valid 
approximately. It is found that the form factors satisfy a relation consistent with the prediction of 
soft collinear effective theory (SCET). In the effective range we investigate behaviors of the form 
factors and differential decay widthes and compare our calculations with the observations from 
other approaches. 
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I. INTRODUCTION 



With numerous scalar meson states being discovered experimentally, most of efforts have 
been devoted to a study about their inner structure and how they are classified. How- 
ever, much controversy persists regarding the underlying components of them. Currently 
one of our main concerns is that these scalar particles can whether or not be described 
consistently in a quark picture. Recently, from a survey of the accumulated experimen- 
tal data the two possible scenarios are suggested , where the scalar mesons below 
and above 1 GeV are assumed to enter their respective nonets in two different ways. In 
scenario 1, there are the two scalar nonets formed by the two quark bound states. One 
contains, as the lowest lying scalar states, the isoscalars cr (600) and / (980), isodoublets 
(k + (800),k°(800)) and (k°(800), «r(800)) and isovector (4(980), ag(980), % (980)). The 
other is made up of the corresponding first excited states: the isoscalars /(1370) and 
/o(1500), isodoublets (K* + (U30) , K*° (U30)) and (K*° (U30) , K*~ (U30)) and isovector 
(do (1450), aQ(1450), <2q (1450)). In scenario 2, those scalar states below 1 GeV are taken 
to be the members of a four-quark nonet, while /(1370), /o(1500), a (1450) and Kq(1430) 
are treated as the lowest lying two-quark resonances and arranged into another nonet, with 
the corresponding first excited states between 2.0 ~ 2.3 GeV. 

Although now we are not able to discriminate among all the existing schemes for the scalar 
mesons, the above two are intriguing in that they can provide us with a ground to make a 
systematic study on the scalar mesons. In such assignment scenarios, an investigation has 
been made into the related decay constants and light-cone distribution amplitudes (DA's)[lj]. 
More importantly, to gain insight into the scalar mesons some of the B decays involving them 
have been explored in the same context. In Ref. jl, 2|, the hadronic decays with a scalar final 
state are discussed in detail in the framework of QCD factorization, important implications 
being drawn for the properties of the scalar particles. More attentions are paid to the 
semileptonic decays with a potential interest — > SWi, Sll. Especially, one shows a great 
interest in the knowledge of their differential rates, since it is critical, as confronted with 
the coming experimental observations, for acquiring valuable information on ingredients of 
the scalar particles. Unfortunately, among the existing approaches no one can afford the 
task to understand the underlying form factors in the whole regions of q 2 , with q being the 
momentum transfers. An effective range of q 2 , in which the calculations are believable, has 
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even not been specified in literature, the computations being carried out in just a small or 
intermediate kinematical region arbitrarily selected. So the results are less persuasive. 

Superior to the three-point QCD sum rules in evaluating heavy to light meson transitions, 
the LCSR approach, which starts with a two-point correlation function, adopts the operator 
product expansion (OPE) near the light cone x 2 = in terms of nonlocal operators, whose 
matrix elements are parameterized as the hadronic DA's of increasing twist. Such that the 
resulting LCSR for form factors, in addition to having an estimable effective region of q 2 , 
can embody as many long-distance effects as possible involved in the decaying processes. 
However, a better understanding of these DA's is critical to have the calculation more 
reliable. Together with the leading twist-2 DA, in general, the twist-3 ones enter and play 
an important role in a LCSR calculation on the form factors. In the case of the scalar 
mesons, the probe into the twist-2 and -3 DA's has been conducted in the framework of 
QCD sum rules and a DA model, in an expansion form in the Gegenbauer polynomials, has 
been formulated, but with a sizable error in some of the model parameters. To try our best 
to reduce uncertainty in LCSR calculation from the long distance parameters, a practical 
improvement scenario has been worked out with its validity examined and confirmed, in 
which a chiral correlator is so chosen that the twist-3 DA's make no contribution [3j. In the 
present work, we intend to apply the same trick to revaluate the semileptonic transitions 
Br s \ Sl&i, Sll, in the two quark picture for the scalar mesons. We will work in the effective 
regions required by the OPE validity and with the two different scenarios aforementioned, 
and calculation is to be performed at leading order in a a . 

The paper is organized as follows: In the following section, we present the correlation 
functions with a chiral current and use them to derive the LCSR for the form factors for the 
— >• S transitions. The discussion and comment are made on the important inputs-the 
DA's and decay constants of the scalar mesons, in Sec. III. Sec. VI is devoted to a detailed 
numerical discussion about the form factors and differential widths for — > Sl&i, Sll, 
including a numerical comparison with the estimates of some other approaches. The final 
section is reserved for a summary. 
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II. THE LCSR FOR THE B 



S FORM FACTORS 



In the standard model (SM), the semileptonic decays Bu^ — > Sli>i, Sll are induced by the 
following effective Hamiltonian: 

Gf- 

w 

G F aV t * b V ts 



eff 



+ 



V2tt 

+ C 10 s7 M (I - j 5 )bi^j 5 £ 
2m b Cf(m b ) _ 



8i<T IJV q v (l+<y 5 )b£'f£ 



(1) 



Here are the Cabibbo-Kobayashi-Maskawa (CKM ) matrix elements, and C\ eS ^ the Wilson 
coefficients, among which C| ff and C\q are scale-independent for the corresponding operators 
have a vanishing anomalous dimension. C^ s and C| ff are expressed as 

Cf(fi) = C 7 (fi) + C b ^ SJ (ri, (2) 
Cf = CM + F pert ( S ') + Y LD (s'), (3) 

where C b ^. sl (^i) stems from the absorptive part of b — > sec — > 57 rescattering which will be 
neglected here, F pe rt and Yld stand for, respectively, the short-and long-distance contribu- 
tions from the four quark operators [4j, with 



^pert 



h(z, s')C - ~h{l, s')(4C 3 + 4C 4 + 3C 5 + C 6 ) 
~h(0, s')(C 3 + 3C 4 ) + ^(3C 3 + C 4 + 3C 5 + C 6 ) 



(4) 



C = 3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 , and 
s') 

h(0,s') 



4 2 I l n 

--ln^ + ^— + -x (2 + x)|l - x\ 1/2 

9 27 9 9 V Jl 1 



2 arctan ■ 



m for x = Az 2 / s' < 1 
for x = 42 2 /s' > 1 



m 6 4 4. 

In Ins H — ^7^ 

9 9 



27 9 /i 



(5) 



TABLE I: The values of Wilson coefficients Ci(m b ) in the leading logarithmic approximation in 
Standard Model, with m w = 80.4GeV, m t = 173.8GeV, m b = 4.8GeV[5|. 



Ci 


C2 C3 


C4 C5 


c 6 


C7 Cg 


C10 


1.119 


-0.270 0.013 


-0.027 0.009 


-0.033 


-0.322 4.344 


-4.669 
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where z = m c /mb and s' = q 2 /m 2 . The Wilson coefficients Ci(rrib), listed in TabJU are given 
in the leading logarithmic accuracy. 

Aiming at an evaluation of the semileptonic decays Bt a \ — > SWi, Sll, we need to confront 
the hadronic matrix elements (S(p) \q 2 ^^i b b\B {s) (p + <?)) and (S(p)\q 2 a flu 'f 5 q' y b\B^(p + q)). 
They can be parameterized, in terms of the form factors f + (q 2 ), f-{q 2 ) and fr{q 2 ), as 



(S(p)\q 2W b\B(p + q)) = -2ip,f + {q 2 ) - i[f + (q 2 ) + f-{q 2 )]q^ (6) 

■ -hiq 2 ) 
m B + m s ' 



(S(p)\q 2 a^ 5 q v b\B(p + q)) = [2p,q 2 - 2q,(q ■ p)] ^ , (7) 



where the B^ s \ mesons are signified by B for short. The relative form factors could be 
calculated in the LCSR. Instead of the correlation functions used in Ref.Jfjj], we would like 
to consider the following two correlators, with the T product of chiral current operators 
sandwiched between the vacuum and one on-shell scalar meson state [TJ]: 



n>,g) = • J d 4 xe^(S(p)\T{ q 2 (x)j„(l - 75 )6(x), 6(0)*(1 - 7s)«i(0)}|0>, (8) 
5 M (p,g) = i [ d 4 xe l(}X (S(p)\T{ g 2 (a;)<v(l + J 5 )q u b(x), 6(0)i(l - 75 ) 9l (0)}|0), (9) 



where q±,q 2 denotes the light quark field. 

The hadronic representations for them are easy to achieve, by inserting between the 
currents a complete set of resonance states with the same quantum numbers as the operator 
6(0)2(1—75)^1(0). On the desired pole contributions due to the lowest pseudoscalar 5-meson 
are insolated, we obtain the hadronic representations: 

jr h( s (S(p)\q2^1 5 b\B(p + q))(B(p + q)\h l5qi \0) 

m 2 B -{p + qY 

| {S{p)\q 2l ,{\ - lb )b\B h {p + q))(B h (p + q)\U(l - 75 ) gi |0) 
^ m 2 B -(p + q) 2 ' 

{S(p)\q2a,u{l + l,)q y b\B{p + q)){B(p + q)\U 1 ,q l \Q) 
m 2 B — (p + q) 2 

{ y (S(p)|g 2 <y(l + l 5 )q u b\B h (p + q))(B h (p + q)\bi(l - 75 ) gi |0) 
^ m 2 B - (p + q) 2 

It should be stressed that the correlation functions receive contributions from the scalar 
resonances included in the intermediate states-B fc |7], in addition to the higher pseudoscalar 
ones, and the ground-state scalar meson is a bit lighter than the pseudoscalar resonance 
lying in the first excited state. 



n*(p,g) 



With the definitions of £>-meson decay constant (_B|6i7 5 gi|0) = ™ B l^ rib and Eqs.flB]) and 
((7|), the phenomenological representations of the correlation functions read 

= m% -t + ™T+^ [2/+(g2) ^ + (/+(g2) + ^ ^ 

i r 2p^( g K + (p^( g )+ P ^( g ))^ 

7T.L *-(p + <?) 2 ' {U) 



m z B - [p + q) z m qi + m b m B + ms 



i r. /r^^-^fa-rii (13) 



7T 7 ao s - (p + g) 5 

Here we have replaced the summations in fflUj) and (fTTl) with the dispersion integrations 
starting with the threshold Sq near the squared mass of the lowest scalar S-meson[7|. The 
spectral densities can be approximated as, by invoking the quark-hadron duality ansatz 

pl,-As)=p Q + c _ D T (s)6(s-s ). (14) 

The QCD spectral densities _^p(s) can be derived by calculating the correctors in QCD 
theory. To this end, we work in the large space-like momentum regions (p + q) 2 << m 2 . 
for the bqi channel and a larger recoil region of the decaying 5-meson as given later, which 
correspond to the small light-cone distance x 2 ~ and are required by the validity of the 



OPE 



10J. Considering the effect of the background gluon field, we can write down a full 



6-quark propagator 

(0|T6(z)6(0)|0) (15) 

Here is the gluonic field strength, g s denotes the strong coupling constant and S (x,0) 
expresses a free 6-quark propagator 

^.0)--/^-i±^. (16) 

The large virtuality of the underlying heavy quarks makes it sound to neglect the contribu- 
tions of soft gluon emission from the heavy quarks, which, in fact, is just a twist-4 effect. 



In this accuracy and 
the traditional LCSR 



eading order in a s , we find that as contrasted with the results of 
6|, only the nonlocal matrix element {S(p)\q2{x)^f^qi(0)\0) remains, 



while those concerning the nonlocal operators q 2 (x)qi(0) and qz(x)(j ^qifi) cancel out. As 



usual, applying the light-cone OPE to the matrix element (5'(p)|g2(^)7^5 , i(0)|0), we could 
be led to the leading twist-2 DA's of the scalar mesons $5(14,//) as defined in [lj. We are 
going to return to this point in the following section. Now the light-cone OPE forms for the 
correlators can be written as follows: 

- ^Sl iu <^h? ' (17) 

b? d m = -^-^■^J^ ^^upr- - (18 » 

We would like to convert them into a form of dispersion integration in order to facilitate 
the ensuing subtraction of the effect of the higher resonances and continuum states in the 
phenomenological representations ( Tl2|) and ( fl3l) . To this end, invoking the relation m 2 — 
(q + up) 2 = u(s — (p + q) 2 ) we make a replacement of u with s. Matching both the forms of 



the correlators, subtracting continuum contributions and making Borel transformation 
with respect to the variable (p + q) 2 , 



1 1 
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B M*— 7 — : — x, = T7?e m 3 

m B ~ \Q + P) M z 

^ 1 _ 1 c -^| 7 [m^+M(l-M)p 2 -(l--u)g 2 ] (19) 

m b ~ (l + U P) 2 uM 2 
with M 2 being the Borel parameter and mj the scalar meson mass, we get the sum rules 
for the form factors: 

/ + (g ) = -^-^ / A du—e , (20) 
/-((fl = ^^^/ ,1 A^^ J (21) 
fr(q 2 ) = q - L fr A {rn B + m s ) / du^^e A } (22) 



where 



A 



2m| 



s — w| — g 2 ) 2 + 4(m 2 — g 2 )m| — (s — m| — g 2 



2 

ml 



A = ~uMi M + M(1 ~ uWs ~ (1 " + (23) 
We find, as a by-product, that the form factors in question respect the following LCSR 
relations: 

/ + (, 2 ) = -/_( ? 2 ), (24) 
Mg2) = (ms±psl f+{q!y (25) 



Actually, apart from that the same is observed in the LCSR involving a pseudoscalar 



meson, a simple relation is obtained also for the form factors in the vector meson case [13]. 
All these observations, up to the hard-exchange corrections, are consistent with the results of 
soft collinear effective theory (SCET) 12] . Having these relations at hand, in the numerical 
discussion we will put our focus on the form factor f + (q 2 ). 



III. DECAY CONSTANTS AND DISTRIBUTION AMPLITUDES OF SCALAR 
MESONS 

In this section, we give a brief review and discussion on the decay constants and DA's of 
the related scalar mesons, which are the basic inputs for the LCSR calculation. 

For a light scalar meson in the two quark picture, it could couple to the corresponding 
vector and scalar quark current operators thus we can define its decay constants as (if, 

(£(p)|<fe(0)7 M ?i(0)|0) = Ppfs, (26) 
(S(p)|fc(0)(&(0)|0) = m s f s . (27) 

It is observed readily that the decay constants fs and fs are scale independent and de- 
pendent, respectively. The neutral scalar mesons like a[j and /o( if considered purely a ss 
bound state) cannot couple with a vector current operator owing to the charge conjugation 
invariance or conservation of the vector current and thus we have 

//o = fao = 0. (28) 

For the other scalar mesons, the decay constants fs and fs are connected by equation of 
motion 

fs = fisfs, (29) 

where 

= 7-v (30) 

m 2 {n) - mi(/i) 

the running quark masses rrii(n) respect the renormalization group equation (RGE): 

a s (/V 



rni(fi) = rrnino) ( > ( 31 ) 
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with b = (33 — 2n/)/3, n/ being the number of active quark flavors. The decay constants 
fs hence are either zero or small of order m 2 — m\. 

Similar to the case of pseudoscalar mesons, the twist-2 DA <3>s(w,/z) of the scalar meson 



is defined as [l] 



{S(p)\q 2 (xhMyM = V, \ dw****™*a{u t ii), (32) 

Jo 

with u being the fraction of the light-cone momentum of the scalar meson carried by q 2 and 
u = 1 — u, and obeys the normalization 

du$ s (u,fji) = f s . (33) 

With reference to the DA's of scalar mesons, a few words should be given. From the 
definition of $5(1*,//), the corresponding scalar mesons have to carry a large light-cone 
momentum p + p 3 . Along with the requirement of the OPE validity, such a constrain 
condition demands that we work in a region assigned as, 

< q 2 < (m 6 - m s ) 2 - 2(m b - m s )A Q cD, (34) 

which, to be specific, is < q 2 < HGeV 2 for a scalar meson below 1 GeV and < q 2 < 
8GeV 2 for one above 1 GeV. Also, it is important to realize that the DA's of scalar meson, 
strictly speaking, become meaningful just at a scale /i > ms, since the constituent quark of 
the scalar meson is in essence off-shell and in particular, it is far from its mass shell by the 
virtuality of m| as carrying the total momentum of the scalar meson. Considering the DA's 
at a scale below mj means that we are dealing with the situation that these off-shell modes 
are in part or in full integrated out, however, which is meaningless. 

Based on the conformal symmetry hidden in the QCD Lagrangian, $,$(«, /i) can be ex- 
panded in a series of Gegenbauer polynomials Cm 2 (x) with increasing conformal spin as 



$ 5 (u,/i) = fs(fM)6uu 



5o(/x) + ^fi m (/i)C^ 2 (2 M -i; 

m=l 

where Gegenbauer moments B m (fj,), which are scale dependent, are given as 

BM - J, ^%ln\n Sl Cl '^- l) ^ )iU - (35) 
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The scale evolutions of $s( u , fi) are determined using the following RGE: 

4/6 



/sO) = /sOo) 



a s (fi) 

where the one-loop anomalous dimensions is [jjj 



B m (fi) = B m {fio) [ , (36) 



m+l 



7M =C F ll - - 2 + 4 £ ~ 

y (m+ l)(ra+2) ^ j 

with Cp = 4/3. The conservation of charge parity demands an antisymmetric $s(w,/i) 
under the interchange u -H- 1 — u, namely, $s(tt, //) = —$5(1 — u, fi), for the neutral scalar 
mesons of a qq content. Accordingly, for the scalar mesons Oq and /o we could write down 
their leading twist DA's as 

*s(u,n) = f s ^)QuuY,B 2m+ i(fi)C^ +1 (2u~l). (37) 

m=0 

In the two quark picture, it is concluded that the twist- 2 DA's of all the light scalar mesons 
are antisymmetric under the interchange u -H- 1 — u in the flavor SU(3) limit, thus the odd 
Gegenbauer moments dominate in the DA's, forming a striking contrast to the corresponding 
situations of the pseudoscalar mesons where the leading DA of the pion, for instance, covers 
no odd Gegenbauer moments and so is symmetric. Indeed the zeroth Gegenbauer moment 
Bo, which is equal to fig , vanishes in the SU(3) limit. In the following, we will neglect the 
contributions of the even Gegenbauer moments and take only into account the first two odd 
moments. 

To proceed, we must add that the LCSR for the form factor f+{q 2 ) would have a distinct 
scale dependence, due to the absence of the QCD radiative corrections. In such a case, it 



should be in order that we work at the scale fib = y rn% s ~ m l>> which denotes the typical 
virtuality of the underlying b quark. At this scale, the related parameters can be evaluated 
making use of the RGE (1361) with an initial scale fio > ms- As the initial conditions we 
prefer using the QCD sum rule estimates at fi = 1 GeV which though is a bit inadequate 
for the situation involving the scalar mesons above lGeV. The numerical results for fs(fji) 
and Bi^ are collected in TabJITI and 1111} and the shapes of the DA's in the two scenarios are 
illustrated in FigJTJ 
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TABLE II: Decay constants f s and Gegenbauer moments ^ of the twist- 2 DA's &s at the scales 
fi = 1 GeV [jj and 2.4 GeV (shown in parentheses) in scenario 1. 



State 


/" {GeV) 


Bi 


B 3 


a (980) 


0.365(0.465) 


-0.93± 0.10 (-0.59 ±0.07) 


0.14 ±0.08(0.07 ±0.04) 


a (1450) 


-0.280 (-0.357) 


0.89 0.20(0.56 ±0.14) 


-1.38 0.18 (-0.71 ±0.11) 


/ (980) 


0.370(0.472) 


-0.78 0.08 (-0.49 ±0.06) 


0.02 0.07(0.01 ±0.04) 


/o(1500) 


-0.255 (-0.325) 


0.80 0.40(0.51 ±0.28) 


-1.32 0.14 (-0.68 ±0.08) 


k(800) 


0.340(0.433) 


-0.92 0.11 (-0.58 ±0.08) 


0.15 0.09(0.08 ±0.05) 



#5(1430) -0.300 (-0.382) 0.58 0.07(0.37 ±0.05) -1.20 0.08 (-0.62 ±0.05) 



TABLE III: Decay constants f s and Gegenbauer moments 3 of the twist- 2 DA's <&s a t the scales 
fj, = 1 GeV [jj and 2.4 GeV (shown in parentheses) in scenario 2. 



State 


f (GeV) 


Bi 


B 3 


a (1450) 


0.460(0.586) 


-0.58 0.12 (-0.37 ±0.08) 


-0.49 0.15 (-0.25 ±0.09) 


/o(1500) 


0.490(0.625) 


-0.48 0.11 (-0.30 ±0.08) 


-0.37 0.20 (-0.19 ±0.12) 


K*(U30) 


0.445(0.567) 


-0.57 0.13 (-0.36 ±0.09) 


-0.42 0.22 (-0.216 ±0.13) 



IV. NUMERICAL CALCULATION AND DISCUSSION 



We proceed to do the LCSR calculation in the two scenarios with the scalar mesons 
in the two quark picture. For illustrative purpose it is sufficient to take, 
study, the processes: B° ->■ a,[ (980) /<( 1450) W h B° s -> /t+(800)/K * + (1430) W h B° ->■ 
RP (800) /K*(U30) II and B?->J J (980)//o(1500) 11. 



The followine inputs , ll6| will be taken in the numerical analysis: 
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FIG. 1: Leading twist distribution amplitudes ^5 of the scalar mesons in scenario 1 and scenario 
2 at the scale \i = 2.4 GeV . It can be seen that $5 is antisymmetric under the replacement of 
u -H- 1 — u in the SU(3) limit owing to the conservation of C parity. 

G F = 1.166 x l(T 2 GeV- 2 , \V ub \ = 3.96^ x 1(T 3 , 
\Vtb\ = 0.9991, \Vts\ = 41.61±g:|g x lO" 3 , 

m u (l GeV) = 2.8 MeV, m d (l GeV) = 6.8 MeV, 
m s (l GeV) = 142 MeV, m b = (4.8 ± 0.1) GeV, (38) 
m e>M = MeV, m T = 1776.82 MeV, 

m Bo = 5.279 GeV, m Bs = 5.368 GeV, 

f Bo = (0.19 ± 0.02) GeV, f Ba = (0.23 ± 0.02) GeV. 
In the first place, let us make investigation in the context of scenario 1. The numerical 
discussions of the form factors f + (q 2 ) can proceed in terms of the standard procedure for sum 
rule calculations. The threshold parameters sq, which correspond to the masses mf of the 
lowest scalar mesons [7], need to be estimated in a certain nonperturbative approach. 
Using the QCD sum rule result s| for the binding energy difference between the scalar 
and pseudoscalar B mesons in the heavy quark effective theory, we could give reasonably 
Sq° = Sq 3 — 33 ± 1 GeV 2 , which is smaller than the threshold values in the corresponding 
conventional sum rule calculations, with the experimental values of the pseudoscalar B 
mesons. Also, it is possible to determine the threshold parameters in other approaches, 
among which the scenario suggested in [9'] is more effective. The range of the Borel parameter 
M 2 , which is shared by all the sum rules in question, is determined as 10 GeV 2 < M 2 < 



12 



TABLE IV: Form factors /+ and /_ at zero momentum transfer q 2 = GeV 2 in scenario 1(S1) 
and scenario 2(S2) for semileptonic decays Bu\ — > Sl~&i with light-cone sum rules(LCSR)[6j, sum 
rules(SR)[17[ and perturbative QCD(pQCD)[18| approaches. 



Methods 


B° s if *+ (1430) 


B° a+(1450) 


B° s k+(800) 


B° a+(980) 


/+ 


f- 


u 


f- 


/+ 


f- 


/+ 


f- 


This work(Sl) 


+0.10 


-0.10 


+0.26 


-0.26 


+0.53 


-0.53 


+0.56 


-0.56 


This work(S2) 


+0.44 


-0.44 


+0.53 


-0.53 










SR [17] 


+0.24 
















LCSR(S2) [6] 


+0.42 


-0.34 


+0.52 


-0.44 










pQCD(Sl) [18] 


-0.32 




-0.31 




+0.29 




+0.39 




pQCD(S2) [18] 


+0.56 




+0.68 













B— >a*(980) 






B— >f (1500) 




Scenario 1 






Scenario 1 




. s =32 GeV 2 




0.4 


s D =32 GeV 2 




s =33 GeV 2 






s =33 GeV 2 




s =34 GeV 2 






s =34 GeV 2 






0.3 








O 
II 

CM 










O" 0.2 








0.1 





10 11 12 13 14 15 10 11 12 13 14 15 

M 2 (GeV 2 ) M 2 (GeV 2 ) 



FIG. 2: Dependance of form factors f+{q 2 = 0) for B° a^"(980) and B° -> / (1500) on the Borel 
parameter M 2 in scenario 1 within the LCSR approach at the scale [i = 2.4 GeV. We take the 
threshold sq = 32,33,34 GeV 2 [7] and b quark mass mj, = 4.8 GeV. 

15 GeV 2 . In this interval, the higher states and continuum contribute less than 30% and 
the sum rule results vary by 13 ~ 30% around the central values, depending on the decay 
modes. 

To elucidate our findings for the form factors, we can consider typically the case of the 
B -+ a (980) and B s -+ / (1500) transitions. The LCSR for form factors, f^^ a t^)^ 
and /^°(1500 )(0)) 

are of a good stability against M 2 -varying, as shown in Figj2j For 
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TABLE V: Form factors /+, /_ and for rare decays -> Sll at q 2 = GeV 2 in scenario 1 

in 



(SI) and scenario 2 (S2) with light cone sum rules(LCSR)[ 



quark modelf 



sum rules(SR 



front (CLF) 



FQM)[22], minimal supersymmetric standard model(MSSM) 



20, 



21], light front 



231 ]. covariant light- 



241 ] . covariant quark model(CQM)[25] and perturbative QCD (pQCD) 



lq approaches. 



Methods 


5° -»• if* (1430) 


B° s / (1500) 


B° -»• k°(800) 


B° s -> / (980) 


/+ 


f- 


h 


/+ 


f- 


h 


/+ 


f- 


h 


/+ 


/- 


/t 


This work (SI) 


+0.17 


-0.17 


+0.24 


+0.14 


-0.14 


+0.20 


+0.46 


-0.46 


+0.58 


+0.44 


-0.44 


+0.58 


This work(S2) 


+0.49 


-0.49 


+0.69 


+0.41 


-0.41 


+0.59 














LFQM [22] 


-0.26 


+0.21 


-0.34 




















CLF [24] 


+0.26 
























SR(S2) [20] 


+0.31 


-0.31 


-0.26 




















SR [21] 




















+0.12 


-0.17 


-0.08 


LCSR(S2) [6j 


+0.49 


-0.41 


+0.60 


+0.43 


-0.37 


+0.56 














LCSR [19] 




















+0.19 




+0.23 


pQCD(Sl) [18] 


-0.34 




-0.44 


-0.26 




-0.34 


+0.27 




+0.29 


+0.35 




+0.40 


P QCD(S2)[18] 


+0.60 




+0.78 


+0.60 




+0.82 














CQM [25] 














+0.40 












MSSM [23] 


+0.49 


-0.41 


+0.60 





















simplicity, throughout the numerical investigation we give only the central values of the 



sum rule results, corresponding to M 2 = 12GeV 2 and Sq° 



B0 



33 GeV 2 . Then we 



have the observations /. 



B°-»<4 (980) 



(0) = 0.56 and /. 



BO^/o(1500) 



(0) = 0.14. Furthermore, use 



of the relations ([24]) and ([25]) leads to /. 



B°-+aJ(980) 



(0) = -0.56, / 



BO^/o(1500) 



(0) 



-0.14 and 



f^ Ml500 \0) = 0.20. Within the LCSR allowed kinematical regions, ff^ +a °( 9S0 \ q i) an d 
/+ s_5 "^ o( ' 1500 ' ) (g 2 ) as a function of q 2 are depicted in FigJHJ along with those corresponding 
to the other modes. The behaviors of f B ^ s (q 2 ) and fT^ S {<l 2 ) are understandable likewise 
with the relations (1241) and f[25l) . Additionally, for a complete understanding of dynamical 
behaviors of the B — > S transitions at the largest recoils, one can be referred to Tab JI VI and 
TabJVl where we collect the present LCSR results for the form factors /_^^^(0) in all the 
cases and the predictions of other approaches for comparison. 
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FIG. 3: Dependance of B^ — > S form factors on the transfer momentum q 2 in scenario 1 within 
the LCSR approach with the scale [i = 2.4 GeV, threshold parameter sq = 33 GeV 2 and Borel 
parameter M 2 = 12 GeV 2 . 
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FIG. 4: Differential decay widths of the semileptonic B — > Slui decays as functions of q 2 in scenario 
1. Here I = e, fi in the left diagram. 



15 




FIG. 5: Differential decay widths of the rare Br s \ — > Sll (I = e, /i) decays as functions of q 2 in 
scenario 1. 
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FIG. 6: Dependance of -B( s ) — > S form factors on the transfer momentum q 2 in scenario 2 within 
the LCSR approach with the scale [i = 2.4 GeV, threshold parameter sq = 33 GeV 2 and Borel 
parameter M 2 = 12 GeV 2 . 

It is manifest that there is a sizable numerical difference in the form factors between 
the transitions to the ground states and to the excited ones. To make it clear, we go back 
to the LCSR expressions for the form factors. We observe that the DA's <3>s(w,/i) make 
contribution only in a smaller region of the momentum fraction u ranging approximately 
from 0.8 ~ 1 at q 2 = 0. The light quark from the heavy quark decays prefers transferring 
to the region close to its kinematical end-point to build a bound state with the spectator 
quark of the decaying heavy meson, which is the so-called Feynman mechanism, that is, 
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FIG. 7: Differential decay widths of the semileptonic B — > Slvi decays as functions of q 2 in scenario 
2. Here I = e, [i in the left diagram. 
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FIG. 8: Differential decay widths of the rare Br s \ — > Sll (I = e, fj,) decays as functions of q 2 in 
scenario 2. 

soft-exchanges predominate over hard-ones in the decay process. Referring to FigJU one 
finds that in that subregion the DA's behave quite differently between the scalar objects 
below and above 1 GeV. For the scalar mesons below 1 GeV, in the whole subregion their 
DA's turn out to be negative and hence make a constructive contribution to the sum rules. 
A different situation manifests itself as scalar mesons involved are heavier ones: the DA's 
contribute constructively in one part of the subrange but do destructively in the other. That 
the two effects cancel out to a large degree leads to a form factor in magnitude much smaller 
than those for the ground states. Physically, this indicates that for a given q 2 , as with the 
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former situation the decaying B mesons have a larger energy release in the latter one. 

In the same picture the B — > S transitions have been explored in the several approaches, 
such as the pQCD [18], QCD sum rules 17, 20, 21] and LCSR Q, Q. It is interesting to 



confront our results with some of the previous studies. In what follows, wherever a result of 
any other approach is referred, it should be understood that we have, if necessary and pos- 
sible, converted it into that in the present convention. Application of the LCSR is enforced 
to B decays to a scalar final state by taking the B s — > / (980) semileptonic processes as a 
study case in Ref. 19j. The sum rules for the form factors, with the asymptotic forms used 
for twist-3 DA's, give /^'^'( 980 )(o) = 0.19 and /^ /o{980) (0) = 0.23 subject to an uncer- 
tainty estimate omitted here. Counting QCD next-to-leading corrections, which is estimated 
roughly based on the observation of the LCSR calculation for the B — > n transitions, the 
above results are modified to /^^^(O) = 0.24 and /s s ^./ (0) = 0.31, about 45% less than 
our calculations. The reason for the sizable differences is mainly use of the different inputs 
for the decay constant // ( 980 ). The different scales are taken for the leading and the sub- 
leading twist DA's as important inputs, which would have, of course, an impact on accuracy 
of result. Using the same inputs, the two evaluations are found to be consistent with each 
other and that of QCD sum rules. The pQCD approach predicts, for the decay modes to the 
scalar ground states, that the form factors are a bit smaller in magnitude but within error 
comparable with the present calculations, and have the approximately same value as in the 
case of the first excited states, a result quite other than our predictions. It is not difficult 
to understand for heavy-to-light transitions, because the pQCD approach accords with the 
hard-exchanges mechanism, and the resulting form factors rely on the behaviors of the DA 
of light meson in the whole momentum region accessible for the constitute quarks. 

All the approaches mentioned above are no doubt applicable in the kinematical region 
near the largest recoil for calculation of the from factors. Nevertheless, no decisive region of 
f, in which these approaches work well, has been provided in the existing applications to the 
B — > S transitions. In the LCSR calculation |6|], the form factors are artificially limited to 
the range < q 2 < 15 GeV 2 , which seem somewhat large against our estimate, and then the 
results are fitted to a dipole model for having an understanding of the behaviors of the form 
factors in the whole kinematically accessible region. The same way is adopted in the pQCD 



calculation 



18| to extrapolate the results for the form factors from small q 2 range to large 



one. Although such a extrapolation manner is phenomenologically extensively assumed, 
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caution should be taken when one applies it to the present case. First of all, we have no 
theoretical justification for doing so. The pole models are believed to be suitable merely for 
description of those form factors corresponding to q 2 near the squared pole masses m^ le , 
however, for the present B — )■ S transitions the m 2 ole are far away from their kinematical 
regions. On the other hand, if the work region for an approach can not be assigned effectively, 
choosing different fitting regions would lead to different results. Hence, it is questionable 
to use a pole description to get an all-around understanding of q 2 dependence of the form 
factors for B — >■ S transitions. Taking this into account, we prefer calculating in the effective 
regions rather than in the whole kinematical range. 

Now, we are in a position to look into the differential decay rates for the B — )■ S semilep- 
tonic decays, which are expressed as 



dT 
dq 



2 ^B( s) Slvi) 



G 2 F \V u b\ 2 q^ 



m 



(q 2 



m 



2\2 



1927r 3 m| 



(q 2 ) 2 
(m 2 + 2q 2 )(q 2 



m%-m 2 s - q 2 ) 2 _ 2 

2 m S 



<r v ^9 

(m B - m s ) 2 ){q 2 - {m B + m s ) 2 )fl(q 2 ) 



+ 3m?(m 2 B — m 2 q ) 2 



dT 
dq 



— 2 {B {S) SIT) 



G 2 F \V tb V ts \ 2 m%a 



2 

em 



i _ 

S 



U(q 2 ) + 

1/2 



1536tt 5 

where mi denotes the mass of a final state lepton, and 



q 2 


m\ 


-m 2 s 


2n 

s 


\ 3/2 



f-(q 2 



,1/2. 



(39) 
(40) 



q 2 /m 



Bi 



T[ = mf/m 



Bi 



Vs 
5s 



(1 _ rs f _ 2s (1 + r s ) + s' 
C 7 f T (q 2 



r s = m 2 s /m 2 B , 



+ \C w f + {q 2 )\ 



<%"U (q 2 ) - 2- x 

6 |C 10 | 2 { [2 (1 + rg ) - s] \f + (q 2 )\ 2 + (1 - r 5 ) 2Re[f + (q 2 ) fl(q 2 )] + a \Uq 2 )f }. 



In the respective effective regions mf < q 2 < (mg — ms) 2 and 4mf < q 2 < (m B — ms) 2 , 
we assess the distributions of the differential rates for the Bu\ — > SWi and Br s ) — > Sll, with 



the results displayed in FigJU and FigJSJ where we have set m e 



m,. 



0. In the case of 



Br a ) — > Sll there appears a discontinuity at q 2 = 4m 2 stemming from the function h(z,s'). 
The differential decay rates for the Br s \ — > St + t~ are incalculable in the present approach, 
for the dilepton threshold 4m 2 is beyond our work regions. It is shown that our calculations 
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and the predictions of pQCD [18( are comparable with each other, although they are based 
on two different dynamical schemes. 

As the scenario 2 is adopted, an analogous LCSR analysis can be made in principle, 
however a complete discussion is not practicable at present, due to little knowledge of the 
4-quark scalar states below 1 GeV. Along the same line as above, we can assess the semilep- 
tonic decays of Bt s \ to a scalar above 1 GeV, which is viewed as a two quark ground state. 
The sum rules show the same Borel interval as in the case of scenario 1. The variations of the 
form factors f + {q 2 ) with q 2 are exhibited in Figj6j and at the largest recoil, a summary of 
the numerical results for the form factors involved, including some of the previous estimates, 
is given in Tab JIVI and TabJVl Comparing the sum rule calculations between the scenarios 
1 and 2, we see that in the the latter case the form factors f + (q 2 ) have a central value 
between 0.40 ~ 0.70 in the effective regions, depending on the decay modes, and hence are 
less sensitive to q 2 than in the former case in which there is a large numerical range from 
0.10 ~ 0.60. The present evaluations of f + (0), which show a better agreement with the 

nn 

conventional LCSR calculation |6|, I19j . are a bit smaller than the numerical observation in 



pQCD 18j, and meanwhile are large numerically in comparison with the calculation of QCD 
sum rules in both the B — > K* and B s — > K* situations, especially our result turning out 
to be about twice as large as that of QCD sum rules in the latter case. 

The resulting differential decay rates, as exhibited in FigJT] and FigJHl have a behavior 
other significantly from what is observed in scenario 1, with the remarkably different QCD 
dynamics embedded in the form factors between the two scenarios. Once these scalar mesons 
above 1 GeV are clearly identified to be, purely or mainly, the two quark bound state, 
this result might help to distinguish between both the pictures for them, as the future 
experiments become accessible. In addition, the distribution shapes, which are demonstrated 
by the differential rates for B^ — > SWi in FigJTJ are compatible with the LCSR calculation. 

The decays to the scalar meson below 1 GeV, despite theoretically little accessible for 
the moment, could be discussed qualitatively. In the four quark final states there is a quark- 
antiquark component from the annihilations of emitted gluons in the decaying processes, 
which gets the transitions highly suppressed. Consequently we may deduce that in scenario 
2 the related form factors are of a small numerical value with respect to the results in the 
two quark picture. 

Finally, we should point out that all the above discussions can not be generalized to 
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decays to a scalar meson, because of the fact that the decaying mesons have a recoil energy 
not large enough to make LCSR applicable, in their decaying processes. 

V. SUMMARY 

We have presented a LCSR computation on B^ — > SWi, Sll at leading order in a s , in 
the two quark picture for the scalar mesons with the two different scenarios. A correlation 
function with chiral current operator is chosen such that the resulting LCSR the form factors 
can avoid the pollution with the twist-3 DA's of the scalar mesons. Applicable regions of 
the LCSR approach are discussed and are assigned reasonably as < q 2 < 11 GeV 2 and 
< q 2 < 8 GeV 2 , for the scalar final states below and above 1 GeV, respectively. Also, 
we investigate the properties of the DA's of the scalar mesons, obtaining an observable 
difference from the case of the pseudoscalar mesons. In the effective regions, the form 
factors and differential decay rates are estimated, with the main findings summarized as 
the following: (1) There exist the relations among the the form factors for the B — >■ S 
transitions, which are in accordance with the prediction of SCET. (2) For the decays to a 
scalar ground state, in the case of scenario 1 the form factors at q 2 = show the numerical 
result much larger those for the first excited sates, and as confronted with the corresponding 
observations in scenario 2, the former seem large in magnitude, but the latter are predicted 
to be small. (3) For the semileptonic processes with the scalar final state above 1 GeV 2 , 
the resulting differential decay rates have a significantly different behavior for the different 
scenarios. Some of them might be beneficial to experimentally identify physical natures of 
the scalar mesons. The present results might be improved as the QCD radiative corrections 
are taken into account, and however they are not expected to change too much from the 
LCSR calculation on the B — >■ it transition {26]. 
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